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ABSTRACT 

We calculate the photon energy dependence of the pulsed amplitude of neutron star (NS) surface modes. 
Simple approximations demonstrate that it depends most strongly on the bursting NS surface temperature. This 
result compares well with full integrations that include Doppler shifts from rotation and general relativistic 
corrections to photon propagation. We show that the energy dependence of type I X-ray burst oscillations 
agrees with that of a surface mode, lending further support to the hypothesis that they originate from surface 
waves. The energy dependence of the pulsed emission is rather insensitive to the NS inclination, mass and 
radius, or type of mode, thus hindering constraints on these parameters. We also show that, for this energy- 
amplitude relation, the majority of the signal (relative to the noise) comes in the w 2-25 keV band, so that the 
current burst oscillation searches with the Rossi X-Ray Timing Explorer are close to optimal. The critical test 
of the mode hypothesis for X-ray burst oscillations would be a measurement of the energy dependence of burst 
oscillations from an accreting millisecond pulsar. 

Subject headings: stars: neutron — stars: oscillations — X-rays: bursts — X-rays: stars 



1. INTRODUCTION 

Accreting neutron stars (NSs) often show nearly coher- 
ent modulations during type I X-ray bursts (see reviews of 
Bildst enll 19981 rS trohmaver & Bildsten 2003) called burst os- 
cillations (Mun o et al.ll200U and references therein). The 
high temporal stability of each NS' s characteristic fre quency 
(within 1 part in 10 3 over years, Muno et al. 2002), along 
with burst oscillations seen from t wo accreting millisecond 
pulsars at their spin f requencies l lChakrabartv et al.l 120031 
IStrohmaver et al.l200 3 ) have led many to conclude that burst 
oscillations are a modulation at the NS spin frequency. Never- 
theless, it has been a long standing mystery as to what creates 
the surface asymmetry in the burst tail, long after any hot spots 
from the burst ignition should have spread over the enti re star 
( Bild stenll99"5l lSpitkovskv. Levin. & Ushomirsk\l2002l) . 

A recently developed and promising hypothesis is that 
burst oscillations are surface r-modes (tHevll 12004ft . In 
this picture the oscillations are created by a retograde 
mode with an observed frequency just below the NS spin. 
As the star cools in the burst tail the mode replicates 
the observed rising frequency. Current theor etica l work 
has focused on calcu lating the frequencies tLed l2004t 



Piro & Bildsten 2005b) and flux perturbations (Hevl 2005; 
Lee & Strohmaver 2005) expected for such modes. Unfor- 
tunately, beside s the highly s inusoidal nature of burst os- 
cillations (Mu no. Ozel. & Chakrabartvl 120021) . which is ex- 
pected for modes, there is little dire ct evidence tha t modes are 
the correct explanation. Piro & Bildsten (2005b) addressed 
this issue by calculating how the mode's properties would 
manifest themselves in the burst oscillations. First, they 
showed that higher persistent luminosity NSs should exhibit 
smaller frequency drifts, consistent with current observations. 
Second, they hypothesized that additional modes might be 
present with such large frequency drifts that they are difficult 
to detect. Though these are promising steps, both predictions 



are directly tied to Piro & Bildsten's (2005b) model that the 
burst oscillation s are a surface wave c hanging into a crustal 
interface wave dPiro & Bildstenll20 05a). What is needed is a 
complementary, and more general, argument of how a surface 
mode should exhibit itself, independent of the specific model 
invoked. That is the goal of this present study. 

A key characteristic of burst oscillations is a larger pulsed 
fraction at increasing energies in the range of 2-23 keV 
(Mun o. Ozel. & Chakrabartvll20 03L hereafter MOC03). This 
is distinct from other pulsing NSs. For example, accretion- 
powered millisecond X-ray pulsars decrease (marginally) in 
ampli t ude between 2 —10 keV (ICui. Morgan. & Tit archuk 
1998; Idallowav et aTJ 120021 Poutanen & Gierlihski 2003). 
MOC03 modeled the burst oscillations as a hot spot on the 
NS surface and nicely reproduced the burst oscillation energy 
dependence. Given the promising possibility that the oscilla- 
tion may be a mode, we perform a calculation of the pulsed 
amplitude as a function of energy for a nonradial oscillation. 
Sim ilar studies have recently b een completed bv lHevll ( 120051) 
and Lee & Strohmaver (2005). Our work differs from these 
in that we are primarily interested in comparisons with ob- 
servations, and thus do not perform an exhaustive parameter 
survey. We also include analytic estimates to explain our nu- 
merical calculations. 

A primary difficulty with comparisons between theory and 
observations is that the mode amplitudes are not predicted 
from linear mode calculations. To overcome this complication 
we follow the example of pulsating white dwarf (WD) stud- 
ies ( Kepler et al1 l2000l) and only fit the shape of the energy 
dependence, leaving the overall amplitude as a free parame- 
ter. Unfortunately, bursting NSs have a drawback with respect 
to pulsating WDs in that the limb darkening of their burst- 
ing envelope is largely independent of photon energy above 
« 1 keV ([Madei 1991). This prevents constraining NS prop- 
erties with burst oscillation data because the inclination, NS 
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mass, M, and radius, R, and even the mode's angular eigen- 
function only alter the pulsed amplitude normalization as we 
show in SI2.3I For this reason, full integrations of the pulsed 
NS emission are well replicated by an analytic result that only 
depends on the NS surface color temperature, T c (eq. l8l). 

We compare our result with the observed energy depen- 
dence of burst oscillations, finding agreement for ksT c m 
2—3 keV, as expected for NSs during X-ray bursts. This sug- 
gests that the burst oscillations are due to a nonradial mode, 
independent of the mode's identification. The excitation and 
nonlinear evolution of the mode is of upmost importance if 
we are to infer NS attributes from burst oscillations. 

In ^2]we calculate the energy dependence of a mode's am- 
plitude and compare it with X-ray burst oscillations in |3] In 
|4] we investigate the optimal photon energy ranges for detec- 
tion of burst oscillations. We conclude in 5j5]with a summary 
of our results and note the importance of measuring the en- 
ergy dependence of burst oscillations from an accreting mil- 
lisecond pulsar. 

2. THE ENERGY DEPENDENCE OF A MODE'S AMPLITUDE 

We first describe our procedure for calculating the pulsed 
amplitude of a surface mode as a function of energy. This 
follo ws previous studies of p ulsed NS emission (for example, 
iPoutanen & Gierlihskil200 3). b ut is incl uded to provide con- 
text for our analytic results in ^2.21 In 32.3l we compare the 
analytics with numerical integrations. 

2.1. Equations for Calculating the Pulsed Amplitude 

We use a spherical coordinate system given by (r,9,<fi) for 
the inertial reference frame of the observer, with its origin at 
the center of the star. The observer sits at an angle 9 = 0. 
Nonradial oscillations are set in a spherical coordinate system 
that shares its origin with the observer's coordinates, but is 
rotated by an inclination angle, i. We denote this by (r, 6', </)'), 
with the pulsation axis, which is coincident with the spin axis, 
at 9' = 0. The cartesian coordinates of the two frames are 
related by 



x =xcos; — zsim 

y'=y 

z' = xsin/ + zcos/. 



(la) 
(lb) 
(lc) 



Gravitational light bending causes photons that reach the ob- 
server to be emitted at an angle a > 9, which is given to high 
accuracy for a Schwarzschild metric by 



1 -COS CH = (1 -COS0) ( 1 



(2) 



where r g = 2GM/c 2 is the Schwarzschild radius ( Beloborodovl 
2002]). 

The equations that we derive adopt a number of simplifying 
assumptions. These are justified because they either are negli- 
gible corrections in the context of pulsed emission from a NS, 
or because they do not affect the energy dependence of the 
pulsed amplitude. We ignore the effects of frame dragging, 
which would modify our results by an amount less than the 
current observational errors (Cadeau et al. 2005 show that a 
Schwarzschild plus Doppler treatment as presented here pro- 
vides sufficiently accurate results in comparison to a full rela- 
tivistic calculation). Lorentz boosting can be ignored since it 
is a < 0. 1 % correction for a NS spin of v = 600 Hz, and as well 
we omit relativistic aberration because it only marginally al- 
ters our results. Finally, we also ignore Doppler shifting from 



the wave motions because the transverse velocity of r-modes 
(~ 10 7 cm s" 1 for an order unity pe rturbation, approximated 
from the results o f iPiro & Bildstenl l2"005b ) is much less than 
c. 

Given these simplifications the observed flux at photon en- 
erg y E is related to the intensity f rom the surface, I{E' ,9' ,<//), 
bv (IPoutanen & Gierlihskil2003l) 

F{E)oc J J S 3 I(E',9',(t)')h(E',cosa)cosadn, (3) 

where h(E', cos a) is the limb darkening function, dSl = 
d cos adefr is the angular element, and E' = Ej (<5 yJ\ — r g /R) is 
the photon energy in a frame co-rotating with the NS surface. 
The Doppler factor is given by 



6 = 1/(1 - /3 sin a sin sin/), 



(4) 



where (3 = 2nRv/ (cy/\ — r g /R) is the equatorial velocity 
(with v the spin frequency). In addition, equation should 
have factors due to gravitational redshift, the NS radius, and 
the NS distance, but we omit these since they cancel when we 
take the pulsed fraction. The integration limits are < a < 
7r/2and0<<?(><27r. 

In general, the limb darkening function can depend on pho- 
ton energy. This is crucial for studies of ZZ Ceti stars, which 
have opacities strongly affected by lines, so that the latitudinal 
quantum numbers can be identifi ed by stu dying th e energy de- 
pendence of the pulsed em ission ( Robinson. Kepler. & Nattier 
119821 iKepler et all 12000). In contrast, bursting NSs have a 
surface opacity dominated by electron scattering. For pho- 
tons with energy > 1 keV the limb darkening is largely en- 
ergy independ ent and well-approximated by h = 0.5+0.5 cosa 
( Madei 1 99 1 ), which is the functional form we assume for our 
calculations. 

Finding the perturbed flux requires perturbing each term 
in the integrand of equation (0 and keeping terms of lin- 
ear order. This results in three integrals to evaluate, which 
correspond to changes in intensity, surface area, and surface 
normal. Since for nonradial incompressible modes the trans- 
verse velocity dominates over the radial velocity (V±/V r ~ 
R/H 1, where H is the scale height i n the bursting layer), 
the latter two changes are neglig ible JButa& Smith! Il979[ 
Robinson, Kepler, & Nather 1982). Using just the integral 
which contains the intensity perturbation, Al, the fractional 
amplitude of the mode is then 



A(E) = 



AF(E) 
F(E) 



J J S 3 AI(E',9' 


', co')/z(cosa)cos ad 


I J d3l(E '- d '> 


<fi')h(cos a) cos adfl 



(5) 



We next relate Al to the mode eigenfunction, which is just 
the temperature perturbation. This relation depends on the 
bursting NS spectrum, which is well-characterized as a dilute 
blackbody with a t emperature g i ven by a color temperature 
T r » ( 1 .4 - 1 . 6)T e sf ( Mad e! 199 It iPavlov. Shibanov. & Zavhnl 
[199T1 lMadejlfl997^l . The change in overall normalization 
does not affect the energy dependence, hence we use I(E') = 
Be'(T c ) and perturb this by setting T c — ► T C + AT, keeping 
terms of linear order in AT, 



AT 



Al _ 9 log/ AT _ x' 

~T~ dio g T c lT~ / -u 



(6) 
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where x' = E' jk^T c . Substituting this result into equation (0, 
the fractional amplitude becomes 



A(E) = 



x'e x AT 



(0 1 , (t)')I(E')h(cos a) cos acM 



I(E')h(cos a) cos ad ft 



(7) 



which can be integrated for any angular eigenfunction 

AT(0',<j>')/T c . 

2.2. Analytic Estimates 

Before we calculate equation @ numerically we simplify 
the integrals so that their energy dependence can be studied 
analytically. There only exists an energy dependence in two 
terms: I(E') and the logarithmic derivative found in equation 
ijS). If these terms contained no angular dependence, they 
could be taken outside of the integrals so that the integrals 
become irrelevant for determining the energy dependence of 
A(E). In principle this cannot be done because E 1 contains an 
angular dependence through the Doppler factor, equation @, 
so that the integrals must be performed numerically. 

On the other hand, if the Doppler shifts are negligible, then 
the energy dependence of A(E) is simply 



A(£)oc 



1 



(8) 



This has high and low energy limits that are useful for gaining 
intuition about the expected dependence on energy, 



A(E)(x 



f E/kaTc, 
1 constant, 



E>k B T cy /l-r g /R 
E < k B T c y/l-r g /R 



(9) 



At high energies the amplitude is linear with energy, while at 
low energies the amplitude is approximately constant. This 
argument shows why a mode will naturally show a larger am- 
plitude at larger energy. It is simply a result of perturbing a 
blackbody spectrum, and the fractional change of intensity is 
much stronger in the Wien tail. Similar results are found by 
MOC03 for their hot spot model when the temperature differ- 
ence between hot and cold regions is small. When the temper- 
ature difference they use becomes large, this can lead to de- 
viations away from a linear amplitude-energy relation at high 
energies (xf > 10), perhaps providing an important discrimi- 
nant between the mode and hot spot models. Unfortunately, 
this is outside the currently observed energy range. 

To find the correction to equation l|8} introduced by Doppler 
effects we expand Al to first order in j3' = j3 sin a sin sin i, 
giving 



AI « (AT) 0=O 



1 + /3' 



2xV 



1 



1 



(10) 



When x' <C 1, this changes the amplitude by a factor 1 + /3', 
which contains no energy dependence. When x' » 1 the am- 
plitude changes by a factor 1 + (3'x' ', so that Doppler shifting 
increases the amplitude at larger energy. 

2.3. Comparisons to the Numerical Integrations 

We now compare our analytic results to numerical integra- 
tions of equation Q. This shows that the analytics reproduce 
the amplitude versus energy relation. For the angular pattern 
of the mode, AT{8', <fi')/T c , we use a buoyant r-mode with 
angular quantum numbers I = 2 and m = 1 (as identified in the 
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FIG. 1 . — The energy dependence of the pulsed amplitude, A(E), for both 
the full numerical integration and the analytic result given by equation |8} 
(thick dashed line). The numerical integrations all use a q = 200, l = 2,m=\ 
buoyant r-mode on a M = 1.4 Mq and i? = 10 km NS. The parameters we 
explore are v = 10 Hz, i = 90° (solid line), v = 600 Hz, ; = 90° (long dashed 
line), and v — 600 Hz, i = 25.8° (dotted line). Though the normalization 
can change drastically for different inclinations (for an example, see Fig. 4 
of Hevl 2005) we renormalize all the results to A(E) = 1 at E/k^T c = 0.1 to 
focus on the shap e of the energy dependence. At high energies, A(E) oc E as 
we show in i|2.2l 



slowly rotating limit, iPiro & Bildstenll2004) . a favored mode 
for reprod ucing the frequency evolution of the burst oscilla- 
tions (lHevl2004l200llPiro & BildsterJ2005bT) . This mode's 
latitudinal eigenfunction is parametrized by the spin parame- 
ter q = 2v I f , where / is the mode frequency in a frame co- 
rotating with the star in units of Hz. A faster spinning NS 
(higher q) results in an eigenfunction that is more concen- 
trated near the NS equator due to Coriolis effects. We fix 
M = 1 .4 Mq and k^T c = 3 keV so that we can concentrate on 
whether other attributes of a NS can affect the energy depen- 
dence. The amplitude at a given energy is assessed by calcu- 
lating the time-dependent amplitude and then fitting this with 
a sinusoidal function. 

In Figure^we compare integrations of different spins and 
inclinations, keeping the mode pattern fixed at q = 200 as well 
as fixing M and R. All amplitudes are normalized to A(E) = 1 
at E/lc B T c = 0.1. At low spin, v = 10 Hz (solid line), the an- 
alytic result of equation (|8j (thick dashed line) and the nu- 
merical calculation are practically identical. As the spin is 
increased to v = 600 Hz (long dashed line) the amplitude in- 
crease at high energies, as predicted by equation JlOi . We also 
consider a more face-on orientation for the NS (i = 25.8°, dot- 
ted line), which has the effect of looking like a faster spin NS. 
This somew hat co unterintuitive result has been seen in previ- 
ous studies JHevlll200^ and is due the mode pattern, which 
has a maximum amplitude at latitudes above and below the 
equator. These comparisons show that neither the spin nor 
inclination change the energy dependence dramatically from 
our analytic result. 

In Figure [2] we keep the spin and inclination fixed at v = 



4 



PIRO & BILDSTEN 



(0 
N 

| 10 
o 

0) 
3 



E 

< 

0) 
W 



1 1 1 1 1 1 1 1 1 1 
r /R=n 41 5, q = 200 


| 


r g /R = 0.208, q=200 




r /R = 0.415, q=10 


h - 


_ Analytic Result 




/ : 








// 


// 


// 

/ 


// 
// 

// 
// ' 


/ ocE 


LI I I I I I I I I 


i i i 1 1 1 



0.1 



1 

E/k D T 



10 



FIG. 2. — The same as FigurelTI but with u = 600 Hz and ; = 90° and varying 
q and R for the numerical calculations. Setting R = 20 km (long dashed line) 
has the effect of decreasing the amplitude at high energies (compare this to 
the solid line). Decreasing q from 200 to 10 (dotted line) affects the energy 
dependence very little, showing that it is difficult to constrain the surface 
pattern created by the mode. 



600 Hz and i = 90°, and investigate the effect of changing q 
and R. When we set R = 20 km {long dashed line) the ampli- 
tude of the pulsed fraction decreases at high energies. This is 
because changing R decreases gravitational redshifting so that 
the break between a constant and linearly increasing ampli- 
tude comes at a higher energy (see eq. |9)). We also decrease 
q dramatically (dotted line) but find very little change in the 
energy dependence. This shows that it is difficult to identify 
the mode pattern on the NS surface via the amplitude energy 
dependence. 

3. COMPARISONS WITH OBSERVATIONS 

MOC03 studied the energy dependence of burst oscillation 
amplitudes from 6 different bursting NSs. A total of 51 burst 
oscillation trains were measured, and multiple trains were av- 
eraged to obtain amplitude versus energy relations. Some 
objects were divided into multiple epochs to assure that the 
gain of the Proportional Counter Array (PCA) on the Rossi 
X-Ray Timing Experiment (RXTE) was relatively constant. To 
correctly compare our calculations to their results we must 
weight our pulsed amplitudes by the PCA's effective area, 
A e f[(E), as well as bin the amplitudes over appropriate energy 
ranges. The PCA is composed of five Proportional Counter 
Units (PCUs). Since each of these have approximately the 
same A e g(E), we use that of PCU3 for our integrations. The 
pulsed amplitude is then given by equation (0 with A e ff(E) 
placed within each integrand. Qualitatively, A e g(E) has a 
large wide maximum spread from pa 4- 15 keV with a smaller, 
secondary peak at 34 keV. The binning of the amplitude de- 
pends on the epoch of the observation and is outlined in Table 
2ofMOC03. 

In Figure [3] we compare the calculated amplitudes with the 
measurements of MOC03 (triangles with error bars). For 



the calculated amplitudes we fix M = 1.4 M@, R = 10 km, 
q = 200, and i = 90. The spin is set to the burst oscillation fre- 
quency for that object. This is reasonable since in all current 
mode explanations of burst oscillations the mode moves ret- 
rogra de with respect to th e spin with v 3> / (see discussion in 
iLee & Strohmaverl2b 05). The normalization of the numerical 
calculations are set to maximize the fit for each comparison. 
We consider k^T c = 3 keV (solid lines) as a fiducial tempera- 
ture exhibited near burst peak. Ideally, we should be able to 
constrain k B T c ^l- r g /R by fitting for the break in the ampli- 
tude (providing M/R if T c is known). This is difficult because 
when the photon energy is > k^T c , as is the case for these ob- 
servations, the amplitude is always linear with energy. Nev- 
ertheless, the theoretical calculations show reasonable agree- 
ment with the observations. 

Comparisons to the observations are complicated by the 
fact that the observations are averaged over a range of temper- 
atures throughout the cooling of the burst, so that we should 
consider temperatures in the range of ksT c k 2- 3 keV. If the 
mode amplitude remains relatively constant, then the pulsed 
amplitude in the Wien tail should increase as the star cools 
(see eq.(9))- We test this in Figure|4]for two of the observed 
amplitudes from Figure|3| but in this case calculating the am- 
plitudes for temperatures of IcbT c = 2 keV (dashed lines) and 
3 keV (solid lines). The overall normalization is chosen to 
maximize the fit, but the relative amplitude of the two curves 
in each panel is set by the temperature ratio. The two tem- 
peratures envelop the data, showing that the spread in the data 
may be due to cooling in the burst tail. An interesting fu- 
ture test of our work would be to divide the burst data into 
early and late stages, and see whether the amplitude evolves. 
This expected temperature dependence can be divided out of 
the data to investigate how much the amplitude of the mode 
is changing due to other effects (e.g., changes in the intrinsic 
amplitude, or changes in q). 

4. OPTIMAL ENERGIES FOR DETECTION OF NEUTRON 
STAR MODES 

Since we understand the spectrum of the burst oscillations, 
this can be used to find the optimal photon energy range for 
burst oscillations searches. The pulsed signal is given by the 
total number of pulsed photons integrated over some energy 
range, 

2 



^ ..I 



R \D 



dE 

fobs / A ef f(E) — 



AT 



1 Tc 



I(E')h(cosa)cosadfl, (11) 



where D is the source distance, f b s is the observing time, and 
we have assumed S « 1. The background noise within this 
energy range is estimated from photon counting statistics 

2 



N = 



'1 
R 



f obs / A ei f(E) 



J J I(E')h(cosa)cos adfl 



dE 
~E 

1/2 



(12) 



which is the square-root of the total number of photons de- 
tected, 

To evaluate equations dl 11 1 and dl 2l i we use a blackbody 
spectrum, I(E') = B E i(T c ), with k B T c = 3 keV, and the analytic 
pulsed fraction, ignoring Doppler corrections, since the agree- 
ment is so close between the numerical and analytic results. 
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FIG. 3. — The observed energy dependence of the burst oscillation amplitudes (MOC03, triangles with error bars), in comparison with our full numerical 
calculations using ksT c = 3 keV (solid lines). We label each panel with the NSs name and the epoch of the measurements if that NS was observed over multiple 
epochs. For each calculation the spin is set to that object's burst oscillation frequency and the normalization is set to maximize the fit. 
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FIG. 4. — The observed energy dependence of the burst oscillations am- 
plitudes (MOC03, triangles with error bars), for 2 of the 9 panels from Fig. 
131 These represent a relatively "good" fit (left panel) and a "poor" fit (right 
panel). We compare these with our numerical calculations using k^T c = 2 keV 
(dashed lines) and 3 keV (solid lines). The overall normalization is set to 
maximize the fit with the data, but the relative amplitude of the two curves is 
set by the two temperatures. This shows that the observed spread in that data 
could arise from the cooling of the NS in the burst tail. 



We assume D = 4.3 kpc (using 4U 1728-30 for a fiducial dis- 
tance), f b s = 10 s (one X-ray burst) and AT/T C = 0.025 (which 
replicates the observed pulsed fractions). In Figure[5]we plot 
these S/N calculations using three different forms for A e ff: a 
flat energy response with A e e = 1000 cm 2 (solid line), the A e ff 
of one PCU from RXTE's PCA (dotted line), and the proposed 
A R ff for the Nuclear S pectroscopic Telescope Array (NuSTAR; 
Harrison et al. 2004, dashed line), a future X-ray mission. 
For each we show a series of 3 keV integrations, which are 
then connected with lines to guide the eye. This comparison 
demonstrates that the energy range of « 2-25 keV contributes 
most to S/N. Other than integrating over this range, there is 
little more an observer can do to maximize the opportunity of 
finding burst oscillations. At high energies the pulsed fraction 
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FIG. 5. — The ratio of signal to noise (using eqs. 1111 and II 21 1 expected 
for a NS surface mode during a 10 s X-ray burst. We set k^T c = 3 keV, 
D = A3 kpc, and AT/T C = 0.025, and compare a flat energy response (solid 
line), the A c ff(£) of one PCU on RXTE's PCA (dotted line), and the proposed 
A c ff for the NuSTAR mission (dashed line). Each curve connects a series of 
points, with each point representing an integration over a 3 keV width energy 
bin. 



is considerably higher, but this range is not necessarily better 
since there are so few photons in the Wien tail and because 
A e ff drops at higher energies. Future missions interested in 
burst oscillations searches could mitigate this by having large 
A e ff at higher energies. 

It is especially exciting that NuSTAR may be able to ob- 
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serve burst oscillations. Though NuSTAR's specifications of a 
large A e ft, high energy range (« 5 - 80 keV), and good spec- 
tral resolution (900 eV at 68 keV) are tuned for observations 
of black holes, active galactic nuclei, and supernova rem- 
nants, it also has fast timing (~ 1 ms) which makes it ideal 
for studying burst oscillations. In addition, its high angular 
resolution (w 40 arcsec) coupled with its timing abilities may 
make it useful for identifying accreting millisecond pulsars in 
crowded fields such as at the Galactic center (something be- 
yond RXTE's capabilities). A typical accreting millisecond 
pulsar at the galactic center has a peak flux in outburst 100 
times less than a type I burst, so the persistent pulse could 
easily be found in a ps 1 day long observation. 

5. DISCUSSION AND CONCLUSIONS 

We have studied the energy dependence of NS surface mode 
amplitudes for NS surface temperatures of k^T c = 3 keV and 
compared this with burst oscillations. The observations fol- 
low our calculated trend of a linear amplitude for photon en- 
ergies > kftT c and becoming flatter when the photon energy 
is « ksT c . Unfortunately, there are currently no data for burst 
oscillations below an energy of k^T c . Measuring the ampli- 
tude at such an energy is crucial for fully testing the mode 
explanation of burst oscillations. One must be cautious of in- 
terpretations for E' < 1 keV because limb darkeningbegins 
to depend on energy in this range, so that equation (|8} is no 
longer applicable. However, this also raises the possibility 
for surface pattern identification a t low energies (ana logous 
to what is done for pulsating WDs, Kepler et al. 2000). 

The agreement between our analytic result and the observa- 
tions is promising for explaining burst oscillations as modes, 
but frustrates the ability of using burst oscillations as a tool to 
learn about these NSs. As long as the amplitude of the sur- 
face perturbation caused by the mode is unknown, it will be 
difficult to constrain NS properties. Future theoretical stud- 
ies should work to address this unanswered question. We ex- 
pect all bursting NS modes to show the energy dependence we 
present here, including the oscillation s een in the 4U 1636-53 
superburst (Strohmaver& Markwardt 2002), for which this 
energy dependence was not determined. 

However, other types of NS pulsations need not match this, 
for example millisecond accreting pulsars in their persistent 
emission. This raises the critical question of what is the 
energy dependence of the burst oscillations from these sys- 
tems (SAX J1808.4-3658 [Chakrabarty et al. 2003] and 



XTE J1 814 — 338 [Strohmayer et al. 200311. |pIro~& BildsterJ 
( 2005b) describe a number of differences between the proper- 
ties of burst oscillations from pulsars and nonpulsars. These 
differences may simply be due to deviations in magnetic field 
strength, but they could instead indicate that the pulsar burst 
oscillations are due to a completely different mechanism. 
Measuring their energy dependence would help to settle this 
unresolved issue. 

One topic we have not addressed is the phase lag observed 
for high energy photons in burst oscillations (MOC03). We 
focused on the amplitude relation because of its stronger sta- 
tistical evidence in the observations. For individual oscil- 
lation trains measured by MOC03, only 13 out of 51 ex- 
hibit phases that vary as a function of energy at 90% con- 
fidence. In comparison, 34 of these exhibit some depen- 
dence with energy at 90% confidence (with the remaining 
typically having less counts in their folded profile). Never- 
theless, it is somewhat troubling that our calculations find the 
reverse phase dependence due to Doppler shifts, just as was 
found in previous studies of pulsed em ission from NSs (e.g., 
IWeinberg. Miller. & Lairibll200UlHevill2005l) . This inconsis- 
tency has been cited by many as evidence that a Comptonizing 
corona exists around a NS during an X-ray burst. Without fur- 
ther theoretical studies of the physical limits of such a Comp- 
tonizing corona or further investigations on how robust of a 
property this observed hard phase lag is, it is not presently 
clear how dire this inconsiste ncy is for th e mode explanation 
of burst oscillations. iLee & Strohmaverl (12005) have found 
some parameter space in R and i that exhibit hard lags, which 
may be promising to pursue further. Though their result is for 
a different eigenfunction than we consider here, it does not 
affect our main conclusions since the energy dependence will 
still be independent of the specific mode. 
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